home
***
CD-ROM
|
disk
|
FTP
|
other
***
search
/
Graphics Plus
/
Graphics Plus.iso
/
general
/
visulztn
/
saoimage
/
saoimage.lha
/
irafenv.c
< prev
next >
Wrap
C/C++ Source or Header
|
1990-04-20
|
10KB
|
300 lines
#ifndef lint
static char SccsId[] = "%W% %G%";
#endif
/* Module: irafenv.c (IRAF Environment)
* Purpose: Do things needed to coorinate paths and coordinates with IRAF
* Subroutine: update_wcs() returns: int
* Subroutine: get_fbconfig() returns: int
* Subroutine: set_path_iraf() returns: void
* Unix calls: fopen(), fclose(), getenv()
* Copyright: 1989 Smithsonian Astrophysical Observatory
* You may do anything you like with this file except remove
* this copyright. The Smithsonian Astrophysical Observatory
* makes no representations about the suitability of this
* software for any purpose. It is provided "as is" without
* express or implied warranty.
* Modified: {0} Michael VanHilst initial version 9 July 1989
* {1} Jay Travisano (STScI) VMS changes 10 Nov 1989
* {2} MVH BSDonly strings.h compatability 19 Feb 1990
* {n} <who> -- <does what> -- <when>
*/
#include <ctype.h>
#include <stdio.h> /* define FILE, stderr, fopen(), fclose() */
#ifndef VMS
#ifdef SYSV
#include <string.h> /* strlen, strcat, strcpy, strrchr */
#else
#include <strings.h> /* strlen, etc. for unenlightened BSD's */
#endif
#else
#include <string.h> /* strlen, strcat, strcpy, strrchr */
#endif
#include "hfiles/coord.h"
#include "hfiles/define.h" /* SZ_FNAME, etc. */
#include "hfiles/image.h"
#include "hfiles/imtool.h"
/*
* Subroutine: set_path_iraf
* Purpose: Construct the pathname of a user datafile by IRAF convention.
*/
void set_path_iraf ( filename )
char *filename; /* i: root filename */
{
static char temproot[64];
char *udir, *getenv();
/* were we passed an absolute pathname as input? */
if( *filename == '/' ) {
return;
}
(void)strcpy(temproot, filename);
/* get defined directory in order of preference */
if( (udir = getenv("WCSDIR")) == NULL )
if( (udir = getenv("wcsdir")) == NULL )
if( (udir = getenv("HOME")) == NULL )
/* at this point we should be checking ~user/iraf/uparm */
#ifndef VMS
udir = "/tmp";
#else
udir = "TMP:";
#endif
if( udir[strlen(udir)-1] == '/' )
sprintf(filename, "%s%s", udir, temproot);
else
#if VMS
/* don't need slash */
sprintf(filename, "%s%s", udir, temproot);
#else
sprintf(filename, "%s/%s", udir, temproot);
#endif
}
#ifdef IMTOOL
static char title[SZ_FNAME];
/* static char comment[SZ_FNAME]; */
/*
* Subroutine: update_wcs
* Purpose: Load the screen WCS, if not yet validated, from the user
* wcs file, if any. (wcs = "WORLD COORDINATE SYSTEM")
* UNIX calls: fopen(), fclose()
* Note: File format (two lines):
* image title (imtool header label string)\n
* a b c d tx ty (WCS coordinate transformation matrix)
*/
int update_wcs ( image, coord, frame_number, wcsbuf )
struct imageRec *image;
struct coordRec *coord;
int frame_number;
char *wcsbuf;
{
float xx, yx, xy, yy, xo, yo; /* l: wcs matrix values */
float low, high; /* l: scaling limits */
int w_type; /* l: scaling code */
int tokens; /* l: number of tokens successfully parsed */
FILE *fp;
char label[1024]; /* l: file title */
char wcs_file[SZ_FNAME];
char root_name[64];
int guess_true_file_coords();
void set_path_iraf(), rotate_transform(), set_trans_speed(), show_filename();
void invert_transform(), combine_transform();
if( wcsbuf == NULL ) {
/* get location of WCSFILE from environment */
sprintf(root_name, WCSFILE, frame_number);
(void)strcpy(wcs_file, root_name);
set_path_iraf(wcs_file);
fp = fopen(wcs_file, "r");
/* if environment not correctly set, try Imtool's defaults, then IRAF's */
if( fp == NULL ) {
sprintf(wcs_file, "%s%s\0", "uparm/", root_name);
fp = fopen(wcs_file, "r");
}
#ifndef VMS
if( fp == NULL ) {
sprintf(wcs_file, "%s%s\0", "~/iraf/uparm/", root_name);
fp = fopen(wcs_file, "r");
}
if( fp == NULL ) {
sprintf(wcs_file, "%s%s\0", "/tmp/", root_name);
fp = fopen(wcs_file, "r");
}
#else
if( fp == NULL ) {
char *getenv();
(void)strcpy(wcs_file, getenv("HOME"));
(void)strcat(&wcs_file[strlen(wcs_file)-1], ".IRAF.UPARM]");
(void)strcat(wcs_file, root_name);
fp = fopen(wcs_file, "r");
}
if( fp == NULL ) {
sprintf(wcs_file, "%s%s\0", "TMP:", root_name);
fp = fopen(wcs_file, "r");
}
#endif
if( fp == NULL ) {
(void)fprintf(stderr,
"WARNING: cannot find imtool coordinate matrix - %s\n",
wcs_file);
return( 0 );
}
/* get iraf coordinate transform matrix */
tokens = fscanf(fp, "%[^\n]\n%f%f%f%f%f%f%f%f%d", label,
&xx, &yx, &xy, &yy, &xo, &yo, &low, &high, &w_type);
(void)fclose(fp);
} else {
tokens = sscanf(wcsbuf, "%[^\n]\n%f%f%f%f%f%f%f%f%d", label,
&xx, &yx, &xy, &yy, &xo, &yo, &low, &high, &w_type);
}
if( tokens < 7 ) {
(void)fprintf(stderr, "WARNING: error reading imtool.wcs file\n");
coord->imgtofile.inx_outx = 1.0;
coord->imgtofile.iny_outx = 0.0;
coord->imgtofile.inx_outy = 0.0;
coord->imgtofile.iny_outy = 1.0;
coord->imgtofile.iadd_outx = 0.5;
coord->imgtofile.iadd_outy = 0.5;
coord->imgtofile.add_outx = 0.0;
coord->imgtofile.add_outy = 0.0;
image->fiscaled = 0;
/* update the transform for tracking info */
combine_transform(&coord->buftofile, &coord->buftoimg, &coord->imgtofile);
return( 0 );
} else {
/* copy the name */
(void)strncpy(title, label, SZ_FNAME);
label[SZ_FNAME - 1] = '\0';
image->filename = title;
/* install its coordinate transform matrix */
coord->imgtofile.inx_outx = xx;
coord->imgtofile.iny_outx = yx;
coord->imgtofile.inx_outy = xy;
coord->imgtofile.iny_outy = yy;
coord->imgtofile.iadd_outx = xo;
coord->imgtofile.iadd_outy = yo;
/* fill in the missing add_out params */
if( yx == 0.0 ) {
coord->imgtofile.add_outx = xo - (0.5 * xx);
coord->imgtofile.add_outy = yo - (0.5 * yy);
} else {
coord->imgtofile.add_outx = xo - (0.5 * yx);
coord->imgtofile.add_outy = yo - (0.5 * xy);
}
/* apply any additional user requested rotation */
if( (!image->row_order) || (image->rotate_code != 0) )
rotate_transform(&coord->img, &coord->imgtofile,
!image->row_order, image->rotate_code);
/* compute speedy integer computation parameters */
set_trans_speed(&coord->imgtofile);
/* compute the inverse parameters (complete) */
invert_transform(&coord->filetoimg, &coord->imgtofile, 0.0);
/* update the transform for tracking info */
combine_transform(&coord->buftofile, &coord->buftoimg, &coord->imgtofile);
/* check for subsection in title and compute true coords (if qp) */
coord->imtool_aux = guess_true_file_coords(title);
/* get IRAF/display scaling (uses linear scaling between low and high) */
if( (tokens < 10) || (low == high) || (w_type != 1) ) {
image->fiscaled = 0;
} else {
image->fiscaled = 1;
image->fiscale = (high - low) / (CMS_DATARANGE - 1);
image->fibias = low - image->fiscale;
}
show_filename();
return( 1 );
}
}
/*
* Subroutine: get_fbconfig
* Purpose: Read the IMTOOL startup file to get frame buffer sizes.
* UNIX calls: fopen(), fclose(), getenv()
* File format: configno nframes width height [extra fields]
* e.g., 1 2 512 512
* 2 2 800 800
* 3 1 1024 1024 # comment
*/
int get_fbconfig ( config_number, width, height )
int config_number;
int *width, *height;
{
register char *ip;
register FILE *fp;
int config, nframes;
char lbuf[SZ_LINE+1], *fname, *getenv();
/* 0 chooses the default configuration */
if( config_number == 0 ) {
*width = DEF_FRAME_WIDTH;
*height = DEF_FRAME_HEIGHT;
return( 1 );
}
/* attempt to open the config file. */
fp = NULL;
if( (fname=getenv(FBCONFIG_ENV1)) || (fname=getenv(FBCONFIG_ENV2)) )
fp = fopen(fname, "r");
if( !fp && (fname = getenv("HOME")) ) {
#ifndef VMS
sprintf(lbuf, "%s/%s", fname, FBCONFIG_1);
#else
sprintf(lbuf, "%s%s", fname, FBCONFIG_1);
#endif
fp = fopen(fname = lbuf, "r");
}
if( !fp )
fp = fopen(fname = FBCONFIG_2, "r");
/* if cannot find a config file, return error */
if( !fp ) {
(void)fprintf(stderr, "Error: cannot find imtool configuration table\n");
return( 0 );
}
/* scan the frame buffer configuration file. */
do {
/* read in next line and check for eof */
if( fgets (lbuf, SZ_LINE, fp) == NULL ) {
(void)fprintf(stderr,
"Error: configuration not found in %s table\n", fname);
return( 0 );
}
config = 0;
/* skip over leading spaces and tabs */
for( ip=lbuf; (*ip == ' ') || (*ip == '\t'); ip++ );
/* skip comment lines and blank lines. */
if( *ip == '\n' || *ip == '#' )
continue;
if( !isdigit(*ip) )
continue;
switch( sscanf (ip, "%d%d%d%d", &config, &nframes, width, height) ) {
case 4:
break; /* normal case */
case 3:
*height = *width; /* default to square format */
break;
default:
(void)fprintf(stderr, "imtool: bad config `%s'\n", ip);
config = 0;
continue;
}
/* is this the configuration we want? */
} while( config != config_number );
if( *width < 1 ) *width = 1;
if( *height < 1 ) *height = 1;
(void)fclose(fp);
return( 1 );
}
#endif